The dynamics of quantum phases in a spinor condensate 



m 
o 
o 

(N 



> 

(N 
O 
m 
o 



o 



S. Yi, O. E. Miistecaplioglu, and L. You 
School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430, USA 

(Dated: February 2, 2008) 

We discuss the quantum phases and their diffusion dynamics in a spinor-1 atomic Bose-Einstein 
condensate. For ferromagnetic interactions, we obtain the exact ground state distribution of the 
phases associated with the total atom number (A''), the total magnetization (AA), and the alignment 
(or hypercharge) (Y) of the system. The mean field ground state is stable against fluctuations of 
atom numbers in each of the spin components, and the phases associated with the order parameter for 
each spin components diffuse while dynamically recover the two broken continuous symmetries [U(l) 
and SO(2)] when A'^ and A4 are conserved as in current experiments. We discuss the implications 
to the quantum dynamics due to an external (homogeneous) magnetic field. We also comment on 
the case of a spinor-1 condensate with anti-ferromagnetic interactions. 



I. INTRODUCTION 

Since the observation of Bose-Einstein condensation in 
trapped dilute alkali atoms Q, |^ |^ , the coherence prop- 
erties of the condensate has become the focus of many 
theoretical studies |1,^ _P SHU- As was initially 
pointed out in Refs. the number fluctuations as- 

sociated with the assumption of a coherent state macro- 
scopic condensate wavefunction leads to the "diffusion" 
(or spreading) of the initial phase. For a scalar, one com- 
ponent atomic condensate, this diffusion, the attempt to 
restore the U{1) symmetry of the underline interacting 
system, can be formulated in terms of the dynamics of a 
zero mode, or the Goldstone mode 3- More physically 
meaningful discussions in terms of the relative phases of 
two condensates were studied soon afterwards ^6., _8 , ^] . 
Starting with the remarkable direct observation of the 
first order coherence in an interference experiment |lfl| |. 
direct relationship between the number fluctuation and 
quantum phase dynamics was first observed with a con- 
densate in a periodic potential [Tll |. and more recently, 
in the remarkable Mott insulating state \13t . 

The emergence of spinor condensates (of atoms 

with hyperfine quantum number F = 1) has created 
new opportunities to understand quantum coherence and 
the associated phase dynamics within a three compo- 
nent condensate 0|. In an earher paper we 
have presented important physical insight into quantum 
phase diffusions of a spinor condensate. We have in- 
vestigated the coupled zero mode dynamics as a conse- 
quence of the usual mean field treatment that calls for 
the breaking of two continuous symmetries: a U(l) gauge 
transformation e'^ and SO (3) spin rotations U{a, (3, t) = 
^-iF^a ^-iFyp ^-iF^T -j-j-^g absence of an external mag- 
netic field) [3. 

This paper provides a detailed investigation of the 
quantum phase dynamics for a spinor-1 condensate. As 
before, we will focus on the case of ferromagnetic inter- 
actions, for which the validity of the single spatial mode 
approximation proves to be a convenient starting point 
for the emergence of transparent physical illustrations of 
the zero mode dynamics The fluctuations of atom 
numbers in each of the spin states will be shown to be 



connected directly, in this case, to the fluctuations of 
the pointing direction of the macroscopic spin (or mag- 
netic dipole moment) of the condensate. This paper is 
organized as follows. In section |nl after introducing the 
Hamiltonian for a spinor-1 condensate, we explicitly work 
out the Heisenberg equations for the atomic held opera- 
tors. Following the standard Bogoliubov approach [2l| . 
we obtain the coupled Gross-Pitaevskii equations for the 
mean field of the condensate and the coupled Bogoliubov- 
de Genncs equations for the quantum fluctuations. We 
then generalize the Hermitian operators for the conden- 
sate number and phase fluctuations (of the Goldstone 
mode) as introduced earlier ^] to each of the spin com- 
ponent, and derive their dynamic equations under the 
rotating wave function approximation. In section UTTI we 
perform a detailed study of the stationary phase diffusion 
of a spinor-1 condensate within the single mode approxi- 
mation valid exactly for ferromagnetic interactions. Sec- 
tion IV addresses the influence of an external magnetic 
fleld, and in section V we discuss the same phase dynam- 
ics for a condensate with anti-ferromagnetic interactions. 
We conclude in section IVTI 



II. FORMULATION 



We consider a system of N spin-1 bosonic atoms inter- 
acting via only s-wave scattering 0, [S^, H^, 13 ■ In the 
second-quantized form, the Hamiltonian of our system 
becomes p.9, 22] 



H = 



2y72 



2M 



+ f E / dr4{7^r^]{7^f^k-f,iU^M^, (1) 

where il)j{r) (j = -|-,0,— ) denotes the annihilation op- 
erator for the j-th component of a spinor-1 fleld. The 
external trapping potential Vc-Kt{r) is assumed spin- 
independent as in a far off-resonant optical dipole force 



2 



trap (FORT) so that atomic spin degrees of freedom are 
completely accessible. The pair interaction coefficients 
are Co = 47r;i^(ao + 2a2)/3M andca = 47rn^(a2-ao)/3M, 
with flo (02) the s-wave scattering length for two spin-1 
atoms in the combined symmetric channel of total spin 
(2), and M is the mass of atom. F is spin 1 matrix rep- 
resentation. The interacting terms can be regrouped to 
show that they include respectively self-scattering, cross- 
scattering, and the spin-relaxation p^ . 

From the Hamiltonian one can easily derive the 
Heisenberg equations for the field operators 



' 2M 



.t 



ipo + 2c20oV'+0- 



" 2M 



+ Kxt(r) 



+C2(-0+V'-(- + -0l-0-)-0O, 



2M 



-0- -I- C20+-0O 



+Co(-0+-0+ + Ipli^O + 0l^/'-)V'- 
+C2(-0l-0- + -0^-00 - -0l-0+)0- 



(2) 



Following the Hartree mean-field theory we assume 
there are three 'large' condensate components around 
which we can study the small quantum fluctuations (off- 
condensate excitations) via the Bogoluibov approxima- 
tion 



(3) 



^j(f)^ ViV0,(r)+(50,(r), 



with N the total number of atoms. We note that 
such an approach can become questionable in a multi- 
component system where certain component becomes 
negligibly small and its quantum nature becomes im- 
portant. For ferromagnetic interactions, however, our 
previous studies [23| show that the ground state wave 
functions can always be expressed as cx and hence 
the Bogoluibov prescription can always be applied. By 
substituting Eq. © into we find in the zeroth order, 
the coupled Gross-Pitaevskii equations (GPEs), 

ih^+ [£ + C0/3 + C2(p+ -I- po - P-)] 0+ + C2iV010O, 

1^00 = [£ + C0/9 + C2(/9+ + p-)] 00 + 2c2iV0S0+0-, 
iTl<j)- = [£ + CoP + C2(p- +P0 -p+)]0- +C2A^0+0O, 

(4) 

with C = —ti^V^/2M + Vext{r), the j-th component con- 
densate density Pj{f) = N\(j>j{r)\'^ , and the total conden- 
sate density p(7^ = p+{f)+pQ{f^+p^{f). J dfpj{r) — Nj 
is the number of condensed atoms in the j-th component. 



In the first order, the quantum fluctuations obey the fol- 
lowing set of Bogoluibov-de Gennes equations (BdGEs), 



ih 



^ [£ + co(p + p+) + C2(2p+ + po - P-)] (5^+ 
+Ar(c+0S0+ + 2c2010o)(5Vo 
-Hc_Ar010+5?A_ + c+N(j)\5i}X 

+ C+iV0+0o5?A(^ + iV(c_0+0_ + C20g)^Vl, 
= [£ + Co(p + po) +C2(p+ +P-)](5-0o 

+iV(c+0;0o + 2c20;^0_)<5V+ 
+iV(c+010o + 2c20S0+)<5?A_ 

+c+N(t)+(t>Q5il}X + c+N (1)0(1)^ S^l_ 

+N{co(j)l + 2c2(b+(j)^)d^jl 

= [C + co{p + P-) + C2{2p- + po - p+)]S'>p- 
+N{c+(j)*(j)^ + 2c20;0o)(5Vo 
-t^c_iV0;0_(5?A+ + c+N(j)lS^l 
+c+iV0_0o(5Vj + Af(c_0+0_ + C2(j)l)S4'l, 

(5) 



where we have defined c± = co ± C2 . 

To study the quantum phase dynamics, we introduce 
the number fluctuation operators [j| 



ih- 



dt 



dr 



4)*{f)Sijjj{r} + 0j(f)(5i/'](r^ 



(6) 



Using the GPEs igj) and the BdGEs ® above, we find 
the evolution equations for these operators, 



ih 



dP+ 



ih dPo ^ (9P_ 
= ih- 



2 dt dt 

= C2N j dr{(t>*^(t>l5i>l + 4>X(t>lH- 
+20!;010o(5-0o - h.c). 



(7) 



Thus both Ptot = P++P0+P- and — are constants 
of motion, which are in fact obvious since the Hamilto- 
nian Q conserves the total number of atoms and the 
total magnetization. 

We define the conjugate phase operator Qtot = Q+ + 
Qa + according to 



~ ih I d 



(8) 



The canonical quantization condition [QtotjPtot] — is 
satisfied if the constraint J+-I- Jo-t- J_ = 1 is enforced with 
Jj = J dr\9*{r)(j)j{r) -|- 6'j (f)0*(f)] . From the defining 
equation for Qtot, 



dQt 



dt 



NiiPt, 



(9) 



the Goldstone mode inertial parameter u can be deter- 
mined. We obtain the dynamic equations for the phase 
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functions as, 

ih0+ = [C + co{p + p+) + C2{2p+ + po - p^)]9+ 

+N{c^(t>+(f)- + C2(f)l)e*_ - Nu4)+, 

ifiOa = [C + co{p + po) + C2{p+ + P-)]So 
+N{c+(l)*^<j>o + 2c2(l)*o<l)-)0+ 
+N{c+(l)*_(t)o + 2c2<l>o(l)+)0- 



+c+N(j)+cf)oeX + c+N(l)Q(j)_e*_ 
+7V(co0^ + 2c2(t)+(t>-)e* - Nu(t>o, 
iM- = [C + co{p + p-) + C2{2p^+po-p+)]0^ 

+N{c+(j)*(j)- + 2c20;0o)6'o + c^N<j)l<p_e+ 

+c+N(t)'i0*_ + c+N4>^(l)oO* 

+N{c^(j)+(j)^ + C2(t)l)6X - Nu(j). . (10) 

Combing Eqs. ijSl and pO|l . the dynamic equations of Qj 
can be easily obtained to be 



Q+ = NuP+ + N J df{~ [(c+(/.;0o + 2c2(/.50-) 6* + c+4>X4>lo^ + {c-4>X4>t + C24>1^) S- + c^^*+(t>-e*_] s^p+ 

+ {c+4>l4>+K + c+<t>X4>W+ + 2c24>-<t>oOX) + {c-(t>\<t>*^e+ + c^4>-(t>+o*+ + C24>fe+) s^^ + h.c] , 

Qo = NuPo + N J dr{- [(c+0*?!)+ + 2c20l0o) 6*; + c+(j)X(j)-^e+ + {c+<j)*(j)- + 2c2(j}*+(po) 01 + c+(t)*(j)*_e^] S^o 

+ {c+(f)X(l)oei + c+(j)*^(l>*0o + 2c2<l)*o<l>-e*) 5i^+ + {c+(l)l<j>o0o + c+<l>o<t>-Oo + 2c2<l>*o<l>+0*o) S^j_ + h.c} , 
Q_ = NuP^ + N j dr{- [{c+(t)*_(t>Q + 2c2(/)S0+) 0*^ + c+(l)*_(t)l0a + (c_0;01 + C2(l)f) 0+ + c-(j)*_(t)+0X] SiIj- 

+ {c+4>l(t>-0*_ + c+(t)*_(t)l0- + 2c20;0o6'* ) + {c^<t)*+((>10- + c_0;0_r + C2(t>f0-) 54,+ + h.c] . (11) 
We note that [Qj,Pk\ — ihdjkJj, [Sipji'f%Pk\ — 6jk4'j{r), and [5iljj{r),Qk] ~ —ihSjk0j{r). This prompts us to make 



the rotating wave approximation (RWA) [7|,|9j by assuming that Sipjir) — J2j ^ji^Pj + X]j ^ji^Qji which leads to 
aj{r) — 0j{r)/ Jj and bj{'r) = (t)j{r)/ihJj. Therefore, under the RWA, we obtain Sipj{f) — 0j{f)Pj/ Jj + (pj {f)Q j / ihJj . 
With this result, the complete dynamic equations for the number and phase fluctuation operators become 



Q+ 



Q- 



_Pq _ 
2 

2c2N 

h 

. 2c2N 
h^ 

NiiP, - 



P_ 

i^-o+m" 



jcf,o0o\ci>+4>-y' ^ 

-r+ — z ro + 

Jo 



2\ii 



J_ 



-P- 



2\i ( Q+ 2Qo , Q- 



2N 



[2c_/_+ + 2C+/+0 + 202(00^01'^-^+)' + C2{(j>l\0-0+)'] P+ 



AN 2N 

-— [C+/+0 + C2{(l),M(ii-0+)'] Po + ^ [2C-I-+ + C2{4>l\0+0-)'] P- 



2C2N 
hJ+ 

Qo = NuPo- 



[2(0+0.100^0)" + (0+^-100 



2\rn 



4c2iV 



i(b-0+\^t>rQo 



47V 



C+/+0 + c+I-o + C2(0-0+|0o0o)' + C2{(t>+0-\^oOoy] Po 



4N AN 

+ — [C+I+O + C2(0O^o|0-e+)'] P+ + -r [c+^-O + €2(00^010+^-)'] P- 
J-l_ J— 



+■ 



4c2N 



[>OOo\(f'+<1 



2N 



)"Qh 



4c27Vj 

hJo 



(0^i0+0_)"]go + ^(0o0oi0+0-)"g- 



NuP- - ^ [2c_/_+ + 2c+/_o + 2c2(0o0o|0+e-)' + C2(0g|0+0-)'] P- 

[c+/_o + C2{cf>o0o\^+O-y] Po + ^ [2C-I-+ + C2{cbl\0-0+y] P+ 
Jo J+ 



-^(0g|0_0+)"Q+ 



^'^'^-(0+0_|02)"go + ^ [2(0_0+ 100^0)" + (0-^+l'^g)"] Q- 



hJo 



(12) 
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where we have introduced the following shorthand no- 
tation; {fg\hy) = J drf*{r)g*ir)h{f)y{r}, (.)' = Re(.), 
(.)" = Im(.), and = / dr{(t>^e*J{<t>ie*)'. Equation 
p2|l thus completely determines the zero mode conden- 
sate fluctuations. Before finding the corresponding com- 
pact zero mode Hamiltonian in quadratic forms of Pj and 
Qj, we will first attempt to simplify in the next section 
using general properties of 4>jif) and 9j{r) for a ferro- 
magnetic spinor-1 condensate |2(]|. 



III. PHASE DIFFUSIONS AND THE ZERO 
MODES 

In general, the ground state wave functions of a spinor 
condensate can be expressed as 



(13) 



with a common chemical potential ^ for all spin com- 
ponents when the external magnetic field is negligible. 
Based on the minimization for the total energy Eq. 
within the mean field theory, it was shown in Refs. 
|25l l26j | that there exists an important relation among 
ttj, given by a+ + a_ — 2ao — 0. We further proved in 
Ref. '2^ that the ground state wave function for each 
of the spin component takes the same spatial shape for 
(ferromagnetic interactions), i.e. 



(14) 



with the real-valued mode function (j){f) (normalized to 
unity) governed by an equivalent scalar condensate GPE 



[/: + c+iV02(f)] 0(f) = Mf), 



(15) 



of a scattering length 02 (note c+ oc 02). Atoms thus only 
collide in the symmetric total spin F — 2 channel in a 
ferromagnetic state in order to maintain their individual 
spins parallel, nj — Nj/N is the ratio of the number of 
atoms in the i-th spin component to the total number 
of atoms. For any given magnetization A4, rij is given 
explicitly by n+ = (1 ± m)^/4 and no = (1 — m^)/2 with 
m = M/N ^20^]. 

Similarly, we look for phase functions of the form 



(16) 



Substituting Eq. ^ into Eq. ((TUJl, we find that 9j{r,t) 
also evolve in time as To focus on the steady 

state quantum fluctuation properties, we may therefore 
neglect the time-dependent part in (f>j (f , t) and 6j (f , t) 
from now on. The steady state equations for phase func- 
tions now become 

+ co(p + 2p+) + C2(3p+ + po - P-)] S.+ 
+2iV(c+0+ + C2(j}-)(j)Qia + N{2c^(j)+(j)^ + C2(f^o)£.- 

= Nu(j)+ , 

[£ - M + co(p + 2po) + C2(p+ + 2N(f>+(t)^ + p_)] 



+2iV(c+(/)+ + C20_)0oC+ + 2iV(c+0_ + C2<I>+)M- 

= Nu4)Q, 

[C- 11 + coip + 2p_) + C2(3p_ + Po - p+)] ^- 
+2iV(c+(/)_ + C2^+)Mo + N{2c-q}+cj)^ + C20g)e+ 

= iVu0_, (17) 

where S.j{r) — 9j(f)e~^°'^ are now real- valued functions. 
One can easily check that similar to the single mode ap- 
proximation for (j)j{r), can also be expressed as 

^j{r) = rij 9{r) with 6{r) governed by the following 
equation 



[C-i^L + 3c+N<j)^{f)] 9{r) = Nu(l){f), 



(18) 



and u is determined through the normalization constraint 
^ dr(j){r)9(f) = 1/2. Therefore, for ferromagnetic inter- 
actions, the phase functions can be generally expressed 
as 



9,{r) = ^,9{r)e^'^K 



(19) 



We note that Eq. (|18|l is in fact identical to the equa- 
tion satisfied by the phase function of a scalar condensate 
a point easily understood for ferromagnetic interac- 
tions when all spins align along the same direction. If 
we take this direction as the quantization axis, then the 
spinor condensate is essentially a scalar condensate since 
its total magnetization M is conserved. All overlap inte- 
grals in Eq. H12|l are presented in Appendix^ They are 
all real quantities with the choice of phase parameters. 
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FIG. 1: The A'^-dependence of u (in units of Tlu>±) for a *''Rb 
condensate in a cylindrically symmetric trap with ujx = u;y = 
Lj± — (27r)100 (Hz) and ujz ~ \lo±. We have used ao = 
101. 8as and a2 = 100. 4as (as Bohr radius) |2^. 

We now attempt to express the zero mode dynamics 
in terms of its associated Hamiltonian. Since Pj and Qj 
are not canonically conjugated variables, we define p'j = 
Pj and q'j = Qj/ Jj such that they become canonically 
conjugates. The zero mode Hamiltonian that governs 
the dynamics of Eq. H12I) is 



H-7.. 



N 



2n+Od 



2c+ 



1 — 



2c2—]p'^O^e 



+ 



2noOA 



2c+^-^-2c2]p'o'0^ 
no ' 



5 



+4 



+4 



-2c+- 



1 — n_ 



C+ + C2 



n 

1/2- 
1/2- 



- 2c2— I p: 

71- 



-4cop+_p'_0^ 



2r 



(20) 



where p'^ = {p'+,p'^,p'_), q'^ = {q'+,q'Q,q'_), A and S are 
two Hermitian matrices whose elements are easily identi- 
fied from above equation. We also present, in Appendix 
IbI the zero mode Hamiltonian in the form of overlap 
integrals. To simplify notations, we have also defined 
O^g = Jdf<j)'^{f)e^{r) and O^^ = J df(f>^{f). We find 
that matrix B is positive definite since C2 < for ferro- 
magnetic spinor condensates. To show the positive defi- 
nitcncss of A, we decompose u as 



(21) 



where 
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FIG. 2: The A''-dependence of e (in units of huj±) for the same 
parameters as in Fig. Q 



{pn,Pm,Py) and q'^ 



{qN,qM,qY)- Not sur- 



prisingly, this transformation recovers the underline sym- 
metry of a spinor- 1 condensate as discovered recently in 
Ref. [2^. It turns out that p^, Pm, and py represent re- 
spectively the fiuctuations of the total number of atoms, 
the magnetization, and the hypercharge. A is not si- 
multaneously diagonalized by the U, yet the zero mode 
Hamiltonian takes a much simpler form 



e = — y dr9{r) - M + c+N(f>'^{r)] Oir). (22) 

We note that £ = for noninteracting atoms [(t>{r) — 
0(r)]. However, it's nonzero in general and we can prove 
that £ > since 9{r) is the solution of Eq. H18|) . By using 
Eq. (PT|l . we then obtain 



pa 



n+ 



no 



where 




A'^\ Co + ^ 



Co 

4c2m 

^0 - 



It is then simply to verify that A' is semi-positive definite, 
which guarantees that A is positive definite. Thus the 
initial ground state, the mean field symmetry breaking 
state with a coherent condensate amplitude \/N tf>j (r) is 
stable. The associated quantum fluctuations of the atom 
numbers can be studied by the linearization approxima- 
tion Eq. 

We note that B is diagonalized by an orthogonal matrix 



1/V3 1/V3 1/V3 \ 
U = ( l/^/2 -1/V2 , 
1/V6 -2/V6 1/^/6 / 



(23) 



which motivates the introduction of new number and 
phase fluctuation operators according to p"^ = p''^U^ = 



N 



ap% 



bpli + cpy 



-apNPY + PpmPy + IPnPm 



-VQy, 



(24) 



where all the coefficients are given in Appendix[ni By re- 
placing with —ihd/dqj, we obtain an eigenvalue equa- 
tion for the distribution of the various fluctuations in the 
condensate ground state 



E 

'm? 



52, 



M 



-a— p 



dqNqY dqAiqy 
V 2 

r29YV- 



We look for eigenfunctions of the form 
which upon substituting into (|25|l yields 



7 



ay 

dqNQM 



fx 

^dql 



IK 



dqy 



E' 



-^lY ) X = 0, 



(25) 



(26) 



(27) 



with 



E' 



E 



-(5k 



M, 



+ bklj+jkNkM)- (28) 
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FIG. 3: The A^-dependence of fi/A^ for the same parameters 
as in Fig. 



The eigen-equation H27|l is essentially of the harmonic 
type, which can be easily solved and the resulting eigen- 
values of Eq. are 



En 

IT 



i2n+l)h{cr,y/^ +A{kN,kM), 



by the Hamiltonian |^ requires (p^) = {p\j) = 0, which 
leads to (g^) — {qlj) — oo, i.e., the ground state of iJzoro 
has completely diffused phases for and qm', The in- 
herent ferromagnetic interaction, nevertheless, prepares 
a correlated ground state such that both py and qy takes 
a Gaussian form distribution with (pv) = 1/2 A^n and 




FIG. 4: The width Ayp of the ground state distribution for 
qy as a function of no and A'^ for the same parameters as in 
Fig. Eland A = 1. 



where n = 0, 1, 2, . . . and 

A{kN,kM) = h^iakjf + bklj +-/kNkM) 



4c 



We note that A{kN, kM) must be greater than or equal to 
zero due to the positive definiteness of matrix A. There- 
fore for the ground state of Eq. H25II . we must have 
kN = kM = 0, the ground state {n = 0) energy and 
eigenfunction are respectively 



Mqy) = ^-i/*(aAyo)"'/'exp - 



91- 




with 



and 



27V. 



YQ 



1 / £-4c20^ 



1/2 



This ground state of -ffzcro can be easily understood. 
For a system starting with a fixed total number of atoms 
and magnetization, the dynamic conservation of the to- 
tal number of atoms and total magnetization as dictated 



{qy) ~ ?i^AyQ/2, which is in fact the minimal uncer- 
tainty coherent state that is consistent with the symme- 
try breaking mean field assumptions of the ground state 
with three separate macroscopic wave functions for each 
of the spin components. However, it is experimentally 
difhcult to produce a condensate having total fixed num- 
ber and phase fluctuations as specified by such a state. 
We therefore turn to study the number and phase fluc- 
tuations of a general state. We start with the dynamic 
equations for pj and qj , 



PN 
PM 

Py 
qN 
qivi 

QY 



0, 
0, 

-2Nr]qy, 

N{2apN - 
N {jpn 4 
N {apN 4 



- 7PM 4 apy ) 
2bpM 4 l3py) , 
Ppm 4 2cpy) 



(31) 



By noting that both pj\j = pAr(O) and pM = Pa/(0) are 
constants of motion, we can express the solutions as 



x(t) = r(i)x(o). 



(32) 



where x^(i) = [pN{t),PM{t),py{t),qN{t),qM{t),qy{t)], 
X^(0) - bAr(0),PM(0),PY(0),gAr(0),gM(0),gy(0)] is its 

initial value, and the evolution matrix T{t) takes follow- 
ing form 
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1 




-a' {I - cosflt) 



V 



CaNt + hAl^Qa'^ sinQt 
(cNt + hAl^f^a'P' sin Qt 
hAy^a' sinQt 




1 

-(3'{l-cosnt) 
CcNt + hA'^^a' 13' sinm 
IbNt + fiA'^^P''^ sinm 



\ 



cosVLt -sinm/?iA|,o 

hA^^a'siwm 1 -a'(l-cosm) 

?iA|.o/3'sinm 1 -/3'(l-cosm) 

?iAyQsinrit cos lit 



(33) 



The coefficients a', /3', Ca, Cfci ^iid Cc are given in Ap- 
pendix ^ The explicit solutions of the time-dependent 
number and phase fluctuations are given in Appendix IeI 
We note that the general solutions has a simple struc- 
ture: in addition to oscillating terms of the form cos fit 
and sinJlt, the phase fluctuations of qjsi and q^v/ also con- 
tains terms proportional to Nt. This indicates that our 
linearization approximation Eq. is valid only for a 
finite duration. The time-dependent covariance can also 
be obtained to be 

{x,{t)x,{t)) ^Y.'^^k{t)T.ji{t){xk{Q)xi{{))). (34) 

kl 

Finally, let's consider the diffusion of the direction of 
the macroscopic spin /(f) = ^'^■(r)'Fij%pj{r). As we 
have shown in Ref. [23| , independent of the spatial coor- 
dinates, the spins of Bose condensed atoms are parallel 
for ferromagnetic interactions, i.e., it essentially acts as 
a macroscopic magnetic dipole pointing along the same 
direction. 



/o(f) = 7V> 0*(f)F,,0,(f) 



V 1 — cos 9 
— Vl — rri^ sin 
m 



(35) 



where Q = a+ — ao — ctQ — a-. Our initial choice of three 
separate macroscopic condensate wave functions (with 
given phases) for each of the spin components fixes the 
initial direction of the condensate spin direction. As the 
phase spreading dynamics attempts to restore the C/(l) 
symmetries of each of the phases, our initial choice of 
phases will become irrelevant. Under the linear approxi- 
mation ||2Jl the fluctuation becomes 



Sf{t) = / df fir) - fo{f) 



(36) 



or explicitly 

(5/^ = J df{<j)*6ij+ + <j)+S^/jl + (j)*_Si;o + ' 
+h.c.), 

5fy = ^jdr{4>lH+ + ^+Hl + 4>*-Ho + ' 
-h.c), 

6f, = [ df{(j)\5i)+ - ^1(5V- + h.c). 



Ml 



(37) 



Although {5f{t)) = 0, the dynamics of the zero mode 
causes the variance of 5f{t) to be nonzero. It can be 
reexpressed in terms of the (time-dependent) number and 
phase fluctuation operators under the RWA as 



5Ut) = 



Sfyit) = 



COS 9 
\/l — m 
a/I — m'^ sin 9 



VSpNit) — V2mpM{t) 



2h 



V2qM{t) + VGrnqyit) 



— sin 9 
Vl — m'^ cos 



VSpNit) - V2mpM{t) 



V2qM{t) + V6mqY{t) 



2h 

5fS) - ^PM{t), (38) 
where the time-dependent terms are 

V2qM{t) + Vernqvit) ^ \/2gA/(0) + VErnqyiO) 

+ - 5- -VSmpNiO) + V2p 

1 — L J 

We see that 6fz is a constant of motion due to the conser- 
vation of total magnetization. On the other hand, both 
6 fx and Sfy diffuse with time. To calculate the variance 
of Sf{t), we simply take 9 = 0, which yields 

1 



Sf.it) = 

Sfy{t) = 

Sf.it) = 



2h 

V2pAf(0), 



^/3pN iO) - V2mpMiO) 
V2qMit) + VGmqyit) 



(40) 



i.e., with the direction of the initial macroscopic spin con- 
strained in the x-z plane. The results of the variances for 
condensate spin direction fluctuations are listed in Ap- 
pendix IfI 

We can also project the fluctuation of the spin [Eq. 
(|2Hll] onto the spherical coordinates (f, 9, cj)) as illustrated 
in Fig. 13 and noting that 



X sm cos (p + y sm sm (p + z cos 6 
X cos 9 cos (j) + y cos 9 sin (f> ^ z sin ( 
~x sin (f> + y cos <p, 



(41) 



and 



sin 9 — \/\ — rn? ^ 
sin0 = — sin 9, 



cos 9 = m, 
cos (f) — cos 0, 
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FIG. 5: Illustration of the macroscopic condensate spin and 
its fluctuation in spherical coordinates. 



we find 

Sfr{t) = 

Sfeit) = 

suit) = 



\/3pjv(0), 
1 

vl — 



2h 



V2qM (t) + Vernqy (t)] . (42) 



We immediately see that both 5fr{t) = 6fr{0) and 
Sfe{t) = Sfg{0) are fixed due to the conservations of N 
and M. Fm-fhcrmorc, with the use of Eq. H39f) . we find 
that Sffii{t) can be rewritten as 

suit) = SUiO) - ^SfeiO), (43) 



where 



suio) 



m 



2h 



V2qMiO) + \/6mqy(0) 



We then see that Sf^{t) grows linearly with t and the 
diffusion rate is given by 

2eN , 



Rd 



-SfeiO), 



(44) 



which is proportional to 6fg{0). A result easily un- 
derstandable in terms of the single axis twisting of the 
isospin (of the Hamiltonian) of a spior-1 spinor con- 
densate [28l| . 

Given any initial condensate state, its variances in 
atom numbers for each individual spin component and 
their correlations and related phase fluctuations com- 
pletely determine the subsequent fluctuations of Sf^it). 
Pj{0) is in fact nothing but the atom number fluctuation 
of state In general we we note that 



df'ip]{r)ipj{r) = N + \/NP.j+ dfS^p]{r)S^pj{r) 



= n + \/np, 



■^"J p2 



2n\j] 



1 

2' 



(45) 



where s = J drO'^{r) and we have used RWA in deriving 
the last line. One can easily find 



5N] = {N^)'{N,) 

dril?Ar)%pj{r) 



dfip'l{r)tpj(r) 



N{P^) 



(46) 



to lowest order. 

If we use cr^, and try to denote the variances of 
the total atom number (iV), the magnetization (A^), and 
the hyper-chargc iY), we find that 



3cr 



[su{t)Y 

{[Sh{t)f) = ^^{3m'<J% + 2ah). (47) 

To calculate the variance of df^(t), we assume a com- 
pletely uncorrelated (between fluctuations oi N, M, and 
Y) initial distribution in the form of a Gaussian wave 
packet 



ViPN,PM,PY) = (27r) ^/''(CTAfCTA/Cry) 
1 fP% 



X exp 



pIi 



'N 



'M 



(48) 



Since qj ~ ihd/dpj, we find for such a distribution 



(9.(0)g,(0)) 



4 2'5u ' 



{P^{0)PM) - '^S.j, 

{qi{0)p^iO)+P^iO)q^m = 0, 

{q^{0)p,{0)) = 0, iori^j, 
and i,j = N,M,Y. Therefore 



[SUitf) = ([SUiO)f 



(49) 

sfemim 



with 



[sumn = 



1 — 



-^ + 3m^^]. (51) 

If on the other hand, we take a distribution similar to 
Eq. H48|l but with no correlations between the popula- 
tions in the spin components j = 0, — , 

Vij^+M-) = (27r)-3/4(a+aoa-)-i/2 



X exp 



1 fpp 



p'i 



pi 



(52) 



with (t| denoting the corresponding atom number vari- 
ance, we find 



([SU{t)f) - al + ol + a^_ 



[sum 



\n \2 2 I 2 2 

[(1 ~m) a++ m 
+ (l + m)V^l, (53) 
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and Eq. (|Sn|l still holds but with 



'sum 



1 — m'^ r (1 + m)'^ 



16 



(1 - m)2 



4m^ 



(54) 



In the Thomas-Fermi (TF) regime, we obtain most of 
the above results analytically. First we note that the 
ground wave function is 



(c+7V)-1[m - Kxt(r)], if > Kxt(^0, 
0, otherwise. 



,(55) 



where c+ = 47r7i^a2/M and the chemical potential 
S-Who / 15A^a2 \ 



with Who = {i^x^yi^zY^^ and a^o = (^i/Afi^ho)^^^- Sim- 
ilarly, ignoring the kinetic energy term in the Eq. p8|l . 
we obtain 



and 



0, otherwise, 



Vo 



(56) 



(57) 



with 



Vn = 



df ~ 



p-Voxt(r)>0 



3 V 'Iho / 



3/5 



All overlap integrals can be calculated analytically as 



'J' 



/ \ 3/5 



167ra^Q Vl5a2 7 

3/2 



252/5 
14^ 



3/5 

— ] iV-3/5^ 

a2 



The parameter e relates essentially to the kinetic energy 
operator Ef^ = -(T? j^M) J drOV'^e evaluated with re- 
spect to the phase function 9{r). In the TF approxima- 
tion we find 



(9) 



Ne = 2E^ 
Using Eqs. ((T^ and lfTH|) we find 



(58) 



which simplifies to 



2Af 



2y72 



2Af 



IL <\>, (60) 



in the TF limit. Multiplying this Eq. (EOl from left by Q 
or ^ respectively, and integrating over r, we obtain 



{4,6) 



f ) -I- 2c+7VO^e 
2c+N I dr(j)^e 



Nil 
Nil 



(4' 

(4' 



(09) 



(0) 



(61) 



where E, 



-{h'^/2M) fdr(f>\/^(f) is the kinetic 



energy per atom in the condensate, and Ej, — 
-{h^/2M) J dr<tN^e = -{n^/2M) J drOV^cj). Eliminat- 
ing E^j^^^ from above two equations, we then obtain 



E, 



Nil 
M 

Nu 

M 
Nil 

M 



E, 



NiL- 



Nu 
' ~2 

2c. N 



- 2cj 



dr(j)^9 



E. 



-E 
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For a spherical trap, adopting the result of Ref. [2£ 
immediately find 



Ah' 



5MNR^ 



In 



R 



l.Bflho 



(62) 



we 



(63) 



with 



which seems to indicate that Ne cx 7V~^/^, asymptoti- 
cally goes to zero at large N. By all measures, we find e 
to be a very small quantity as compared to 0204,0- Tak- 
ing e as zero in the TF limit, we obtain 



nn = - 



1 / 7 
sV 10 ' 



151/10 



C2 



na 



ho 



3/5 



0-2 



To shed light on the above scaling results based on 
the TF approximation, we have performed extensive nu- 
merical calculations. As an example, we consider a spinor 
condensate of ^^Rb atoms in a spherically symmetric trap 
with a harmonic frequency uj = (27r)100 (Hz) and con- 
sider the regime of iV ~ 10^ — 10^^. It turns out that our 
results indeed confirm the scaling relationship 



and 



Nil 



YO 



2 

5' 

1 rr 
3 V 10' 



(65) 
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in the large N limit. We also find that 

j^W ^ ^-0.3549^ 

in rough agreement with the result of Ref. ; On the 
other hand, we find 



and 



Ne w AE^ 



k ' 



which do not agree with the results above as obtained 
with the TF approximation. A careful analysis as con- 
firmed by our numerical results reveal that although 
/ d'r{Vc^t - Ai + c+N(jP')(lP' is small, / dr{Vc^t - + 
Cj^N(jP')9'^ is not, and approximately we find 



E 



(6) 



(66) 



This leads to the opposite A^-dependence for El. and 



eI'''^; while El"' increases with increasing values of N, 

e\^^ decreases. Thus the identity Eq. H62|l as obtained 
in the TF limit is invalid. 



IV. A NONZERO MAGNETIC FIELD 

Inside a nonzero homogeneous magnetic field {B), the 
atomic quantization axis becomes fixed, along the di- 
rection of i3, or more conveniently denoted as the z- 
axis. Mathematically, it turns out this is equivalent to 
the the requirement of the conservation of magnetization 
M — A^+ — A_. This reduces the S0(3) rotational sym- 
metry to its subgroup SO (2), which is still a continuous 
symmetry. Thus the quantum dynamics of the phase re- 
mains formally the same as discussed before. 

In the linear Zeeman regime, the original Hamiltonian 
(Q) is augmented to 



H' = H - hojiF^ 



(67) 



The Larmor processing frequency is lol = Bfj,B/fi with 
fiB the magnetic dipole moment for state \F = 1, Mp — 
1). In this case, the only change is the phases of the 
spatial wave functions and the phase functions if of the 
following time-dependent form 



(68) 



for j = -l-,0, and — . a+o — aoo = aoo — a_o = still 
holds true for a ferromagnetic spinor-1 condensate. A 
direct consequence of the magnetic field is thus a time- 
dependent relative phase 



e(i) = a+o - aoo + t^L* 



(69) 



When applied to the macroscopic condensate spin, we 
find the spin gains a precession around the z-axis. 

Now the picture of the dynamics of the macroscopic 
condensate spin is clear: besides the precession induced 
by the external B-field, it also diffuses to restore the U(l) 
and SO (2) symmetries of the Hamiltonian. 

The effect of quadratic Zeeman shift can also be simply 
addressed. Apart from an overall shift, it introduces a 
(positive) B-field dependent level shift cx HSbF^ to states 
|±) with respect to state |0). This differential shift causes 
the precessing of the condensate spin to be twisted, i.e. 
the precessing of the |-|-) (|— )) component being slower 
(faster) by (5s- 



V. THE CASE OF ANTI-FERROMAGNETIC 
INTERACTIONS 

Our discussion of the phase diffusion in the previous 
section is based on the proof that the single spatial mode 
approximation is exact for ferromagnetic interactions. 
One might therefore conclude that a similar study for 
anti-ferromagnetic interactions would be difficult as, un- 
like for ferromagnetic interactions, it's difficult to predict 
the ground state condensate structure. In reality, how- 
ever, the quantum phase problem for anti-ferromagnetic 
interactions is much simpler, provided the mean field de- 
scription for the ground state still applies. 

First we note that when Ai 0, the population in 
the 0-th component is zero in the ground state 20, 2^. 
Therefore, the zero mode Hamiltonian H20|) reduces to 



Pi 



N 



zero + / ~ I o T 

Pi f J_ 



ic-P+P- 
J+J- 



(70) 



which is exactly the same zero mode Hamiltonian for 
a binary condensate with no coupling between its two 
components. This has already been studied before 0- 
We note that results obtained in Section ^remain valid. 

Now, let's consider the Ai — case. As was pointed 
out in Ref. |23, the single mode approximation again 
applies in this case, and we can assume 

0,(f) = vH70(f)e-^ (71) 

where = n_ = (1 — no)/2 and total energy of the 
system is independent of the value of riQ. A similar re- 
lation exists among the phases of the three components: 
Q!+ -|- q;_ — 2ao = TT. The spatial profile 0(r) satisfies 
following equation 

[C + coN(f>''{r)] ^(r) ^ fi<f>if), 



(72) 

with a coefficient cq for the nonlinear term. The phase 
function then becomes 



(73) 
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where 6{r) is the solution of the following equation 



we obtain 



[C-fi + 3coiV(/>2(f)] 9{r) = Nu0{r}. 
Following the results of Sect. IIIII we define 



(74) 



N 



dre{r) [C-^l + coAr,/,2(f)] 9{r), (75) 



and the zero mode Hamiltonian becomes 



N 



1 - m)04>e 
e 



2cn 



2C2 

1 - no 



2co p'o'O^ 



2co 



2C2 



(1 - ^0)000 ■ ' ' 1 - no 
4copV-Po'^0e + 4cop'_PoO0e 

2no - 1\ , , 

P+P-Od 



p'lO, 



+4 Co + C2- 

V 1 - "0 

^ C2no(l - na)0^^ 



(76) 



the positive definiteness of matrices A and B can again 
be verified. By applying the same transformation (|23|l . 



N 



ap% 



bp 



M 



cPy 



aPNPY 



(77) 



where all coefficients are now given in Appendix [H] We 
see that fluctuations in N are completely decoupled from 
that of M and Y in this case, a result that also happens 
for ferromagnetic interactions when = 0. The ground 
state of (|76|l is similar to Eq. 1)20(1 except now that 



2N 

2Ny^ = —[C20^^{e + 4c2ngO^ 



11/2 



1/2 



(78) 



The dynamic equations for pi and qi then become 

Pn = 0, 

Pm = 0, 

Py = -2Nr]qY, 

qN ^ N {2apN + apy) , 

qM = 2NbpM, 

qY = N {apN + 2cpy) ■ (79) 

One can easily find the evolution matrix to be 



1 



^a'{l - cosnt) 



CaNt + hA 



YO'- 



■ sin fit 







hAyQa' sin Qt 




1 



2bNt 







cos ilt 









h A^r^a' sin fit 1 





^Ayq sin fit 



1 







- sin flt/hAyg 

-a'il - cosilt) 


cos fit 



(80) 



where 



VI. CONCLUSIONS 



Ca 



a _ (3no - l)e + 8c2nl0^e 

2'c ~ V2(e + ic^nlO^g) ' 

iac-a'^ 3e{e + ic2no0^g) 
= UcoU^e + ^ , o, 



2c 



s + 4c2n^C'^ 



As for 6f{t) we find 



^ / 2no ^6no(l - no) 
i 8 W Pm + sm 6 

1 — no h 



qy, 



. r^ I 2?io , \/6"-o(l - m) 
— smB^ / z Pm + cosfc)-^^ gy, 



1 - no 



= V2p_ 



M- 



(81) 



In conclusion, we have investigated in detail quantum 
phase diffusions of a spinor-1 condensate. When the elas- 
tic atomic interaction is of ferromagnetic type, the struc- 
ture of the ground state is greatly simplified: the spatial 
mode functions of the different spin components are ex- 
actly the same and accordingly we also proved that the 
phase functions are exactly the same. This simplifica- 
tion allows us to construct the zero mode Hamiltonian 
which describes the number and phase fluctuations of a 
condensate while maintaining the conservations of the 
total number of atoms and the total magnetization. We 
have provided analytical results for the number and phase 
fluctuations due to both the quantum phase diffusion dy- 
namics and the initial distribution of atom number and 
phase fluctuations. The structure of these fluctuations is 
rather simple: along with the oscillations due to sin fit 
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and cos Qt terms, quantum phase diffusion terms propor- 
tional to Nt also exist. This suggests that the mean-field 
approach to spinor-1 condensates is only valid for a finite 
duration. Based on these results, we have studied the dif- 
fusion of the direction of a macroscopic condensate spin. 
Our investigation sheds important light on the studies of 
beyond mean field theory quantum correlations among 
Bose condensed atoms. 
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APPENDIX A: OVERLAP INTEGRALS 

In this appendix, we list the overlap integrals as in- 
volved in the dynamic Eq. ifT^ . Using the phase con- 



ventions as introduced in Eqs. 1)14(1 and (|19() . we find 
explicitly 





~ ""i ' 












- {4>l\o+e^) 















APPENDIX B: THE ZERO MODE 
HAMILTONIAN 



The zero mode Hamiltonian is listed below. 



N 



-^u + 2c_/_+ + 2C+/+0 -I- 2c2(<^o0o|</'-0+)' + C2((^^|0+0-)' 

2_Pl 
p2 r 

4P P 4 P P 

+ ° [c+I+o + C2(0o0o|0-^?+)'] + [c+I-0 + C2(0o^o|0+e-)'] 



-^u + C+I+o + c+I-0 + C2(</)_0+|0o^o)' + C2{cf>+9-\(l)o0oy 



+ 2c_/_+ + 2c+/_o + 2c2(0o^ol0+^-)' + C2{(l^l\e+e_y 



J+Ja 
2P+P_ 



J-Jo 

/ ,2| 



(Bl) 



r 



APPENDIX C: COEFFICIENTS IN THE ZERO 
MODE HAMILTONIAN FOR FERROMAGNETIC 
INTERACTIONS 



The coefficients in the zero mode Hamiltonian (|24f) are 



4^0 L 



-(5-t-3m2) 



+204,6 (^3co(l - mY + 3C2(1 - 3to2) 
Zrir, 



/2 

a ^(1 + 3to2)(;._4c2 



2m 



(e - 4C2O00), 



2m 2 r „ / 
7 2\h [£-C20^e{l + 3m^)] , 

Hq V o 

rj = -2,C2nlOif>4,/h^ . 



(CI) 
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APPENDIX D: PARAMETERS IN EVOLUTION 
MATRIX 

The list of parameters in the evolution matrix T(i) 
(ESI are 



46c - /?2 _ 4e 

2c 1 — ' 

2c7 — a/3 2\/6me 



2c 



1 — 



(Dl) 



a 



a 1 + 3m2 



2c 2^/2 ' 

/3 



/3' = ±l = _V3m, 
2c 

4ac-a2 3(l + m2)e 

Ca = ?r:^ = ZTY- + 12(C0 + C2 



2c 



1 — 



APPENDIX E: EXPLICIT FORM OF NUMBER 
AND PHASE FLUCTUATION 



Here we present the explicit form of the solutions for 
the number and phase fluctuations. 



PNit) 
PMit) 

PY{t) 

qN{t) 



qM{t) 

qvit) 
p+{t) 

Po(.t) 
pUt) 

q+{t) 



pn{0), 

Pm{0), 
1 + 3rrf 



2V2 



■(1 - cosnt)pMiO) + %/3m(l - cosm)pM(O) +py(0)cosm - |^|^ 



2 sin rii, 

YO 



Nt I + + 12(co + C2)0^e ) + 
1 — ' 



(l + 3m2)2 . ■ 
—fi^Yo sill 



Pn{0) 



Nt 



2\f%me \/3m(l + 3m^) ^ . n . ^ 

-?iAyQsmS« 



1 — 

2 



2\/2 



Pm(0) 



1 + 3m'^ 1 + 3m^ 

-?i^A|.oPY(0) sin VLt + (7jv(0) ^7^(1 " cos rJi)gy (0), 



2\/2 

2A/6m£ 



2\/2 



V3to(1 + 3to2) . 



z sin f2i 



nA 



I — rn?- 

As 
I — m? 
1 + 3m2 



2\/2 



Piv(O) 



im^fi/^Y^ sinf7i ) pAf (0) - \/3m;i AyoP^ (0) sinf7t + 1771/(0) + Vini{l - cos fit)(jy (0), 



PAr(O) sinili - ?iAyoV3TOpj\/(0) sin Vlt + ?lAyQpy(0) sin fit + qy(0) cos fit, 

sin fit 



3(1 - m^) + (1 + 3™^) cos fit 





4\/3 




3(1 + m2) 


-(1 + 


3m^) cos fit 




2\/3 




3(1 -m2) 


+ (1 + 


3m^) cos fit 



pn{Q) 



1 + TO — m cos fit , , cos fit , , 
-Pm(0) + ^^Py(0) 



V2 



V6 



V6?.A|.„''^(°)' 

J- , ^ . , , 2 cos fit , , 2 sin fit , , 
pjv(O) - \/2to(1 - cosf!t)pM(0) 7^Py{Q) + ^=-— gy(0), 



4V3 



PJv(O) - 



%/6 



V2 
V3?iA2 



V6 



V6?iA2,Q 
sin fit 



VSA^t (L^^ ^ 4(^^ ^ C2)O0e) + (1 + 3to2)(1 + mf sin fit 

\ 1 + m y 
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2^/2e 3to(1-to 2 

— Nt ^ — ^-^HAtr^ sin f7t 

l + m 2V2 ° 



^.2 3(1 - m)2sinf7t 
Pm(0) + ;iA2,o-^ 



-^9w(0) + ^'ZAf(0) 



1 — 6to + 3m^ (1 — m)^ 



2^6 



2V2 



cos fit 



2V6 

qviO), 



PA'(O) 

Py{0) 



q'oit) - 



VSA^t f l^-^e + 4:{co + C2)0^e] + ^hAl^Q{3m'^ - 2m^ - l)sinf7t 



1 — TV? y 

2V2me 3to(1 - to^) . ^ 

-^^ Nt + h Air. ^ — - — ^ sin nt 

1-to2 2\/2 



Pw(0) 



Pm(0) + hAl^o sin»py(0) 
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V3 



1 + Sm^ 3m^ 



2V6 



2V6 



■ cos fit 



\ I — ni 



Pn{0) 



2V2e 3to(1 + to)2 

1-m 2V2 



/X * . 9 3(1 + m)2 sin fit , , 



+ ^9^(0)- ^'Zm(0) + 



1 + 6m + im? (1 + mY 



2V6 



2V2 



cos Vlt 



2V6 

9^(0). 



(El) 



APPENDIX F: THE VARIANCE OF 
MACROSCOPIC CONDENSATE SPIN IN 
CARTESIAN COORDINATES 



In Cartesian coordinates, the variances of the conden- 
sate spin direction is given by 



[5fy{t)f) = 



1 — 

1 — m? 



3(p^(0)) - 2\/6m(pA,(0)pM(0)) + 2rr?{plM) 
^MM + 2N/3m(9M(0)gy(0)) + 3m2(gf,(0))] 



2% 

V2Net 



■[Sm^iplm - 2V6m{pN{0)pMm + 2{pli{0))] 



V2{(qM{0)pM{0)) + (pM(O)gM(O))) - V3m{{qM{0)pN{0)) + (pw(O)gAf(O))) 

n L 

+\/6m(((Zy(0)pjv(0)) + (pAr(O)gy(O))) - 3m^ ((gy(O)pAf(O)) + (pAf(O)gy(O)))] , 
2(pL(0)>. 

I 



APPENDIX G: COEFFICIENTS IN THE ZERO 
MODE HAMILTONIAN FOR 
ANTI-FERROMAGNETIC INTERACTIONS 

The coefficients of zero mode Hamiltonian (|77|l for anti- 
ferromagnetic interactions (with = 0) 



4c2no 



-a 



3no(l-"-o) 3(1 -no) 
2(3no — l)e 16c2rio 



3%/2no(l-no) 3\/2(l - tiq) 
77 = 3c2no(l - no)©^^/?!^. 



0<i 



(1 + 3no)e 



6no(l - no) 



6co 



3(1 -no) 



1 - no 



(Fl) 
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